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Results are presented of an analysis of proton and charged pion azimuthal 

distributions measured with respect to the reaction plane in Au + Au colli- 
sions at a beam momentum of about 11 A GeV/c. The azimuthal anisotropy is 
studied as a function of particle rapidity and transverse momentum for differ- 
ent centralities of the collisions. The triple differential (in rapidity, transverse 
momentum, and azimuthal angle) distributions are reconstructed. A compar- 
ison of the results with a previous analysis of charged particle and transverse 
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energy flow as well as with model predictions are presented. 



I. INTRODUCTION 

Collective phenomena play an important role in heavy-ion collisions, but for a long time it 
was assumed that, at collision energies much greater than 1 GeV/nucleon, only longitudinal 
and azimuthally symmetric transverse radial flow would survive. During the last few years 
the situation has changed qualitatively. Anisotropic, directed as well as elliptic, flow has 
been observed at the BNL AGS and strong indications of elliptic flow at the CERN 

SPS 1^ have been demonstrated. The theoretical understanding of the effect and model 
calculations involving anisotropic flow have progressed significantly; different anisotropic flow 
patterns have been linked to such phenomena as quark-gluon plasma formation, softening 
of the equation of state and mean- field effects during the fireball evolution 



The appropriate tools for flow studies at high energies have been developed [[rii|T2| . It was 



also noticed that anisotropic flow could be important for other measurements, such as two- 
particle correlations |jIB|,p!^. Anisotropic flow has become an essential part of the global 



picture of heavy-ion collisions at ultrarelativistic energies ; it is considered one of the key 
tools to elucidate the dynamics of the collision. 

In the current paper we present results of the analysis of anisotropic transverse collective 
flow of identifled particles, protons and charged pions, detected in the E877 spectrometer 
in Au -|- Au collisions at a beam momentum of 10.8 and llAA GeV/c. The data were 
taken during the 1993 and 1994 AGS heavy-ion runs. Using calorimeter data we reconstruct 
the reaction plane event-by-event, and analyze the particle production with respect to this 
reaction plane. A similar analysis of charged particle multiplicity and transverse energy flow 
presented in displayed a strong directed flow as well as an elliptic flow with the primary 
axis in the reaction plane, i.e. perpendicular to the direction of the "squeeze-out" effect 
observed at lower energies |T^. Here, using the E877 spectrometer data, we apply the same 
procedure to determine the triple differential distributions N / {dy ptdptd<P) of identified 
protons and charged pions. 
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The paper is organized as follows. After a description of the apparatus we discuss 
the procedure of the analysis. In the subsequent sections we present first the results of 
the Fourier analysis of azimuthal particle distributions in different rapidity bins and for 
different centralities of the collision. These results are, to a large extent, independent of the 
uncertainties in the spectrometer efficiency and the detector acceptance. Then the results 
on azimuthal anisotropies and related quantities (such as the mean transverse momentum 
projected into the reaction plane {px)) are discussed. Different flow scenarios are discussed 
from the point of view of the observed dependence of the flow signal on the transverse 
momentum pt of the particle. We compare our measurements with the evaluation of nucleon 
and pion flow derived from the measurements of charged particle and transverse energy 
flow as well as with model predictions (RQMD versions 1.08 and 2.3 |T0| , p!7[] ). 



II. APPARATUS 

The E877 apparatus is shown in Fig. 1. In the E877 setup, charged particles, emit- 
ted in the forward direction and traversing the collimator (—134 < Ohorizontai < 16 mrad, 
— 11 < 9yf,rticai < 11 mrad), are analyzed by a high resolution magnetic spectrometer. The 
spectrometer identifies particles via simultaneous measurements of momentum and veloc- 
ity. The momentum of each particle is measured using two drift chambers (DC2 and DC3, 
position resolution about 300 /im) whose pattern recognition is aided by four multi-wire 
proportional chambers (MWPC). The average momentum resolution is Ap/p ^3% limited 
by multiple scattering. A time-of-fiight hodoscope (TOFU) located behind the tracking 
chambers provides the time-of-flight with a typical resolution of 85 ps. The spectrometer 
acceptance covers mostly the forward rapidity region with transverse momentum coverage 
including pt = 0. Further discussion of the acceptance for different particle species can be 
found in . 

A clean particle identification is particularly crucial for the study of flow since the ex- 
tracted signal may be different for different particle species both in magnitude and in sign. 
A small admixture in the particle sample of other particle species could therefore heavily 
bias the final results. Consequently, in the current analysis, strict requirements are applied 
for particle identification. Using the measured momentum and velocity of the particle we 
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calculate the particle mass. In a plot of momentum vs. the square of the particle mass 
we select entries in the region of ±1.5(Tm2(p) around the peak of the particle under 
study and reject all entries within a 3 a^2(j)) region of another particle species. Here, ijm'^{p) 
represents the standard deviation of the (Gaussian) mass-squared distribution at a given 
particle momentum p. 

The determination of the centrality of the collision and of the reaction plane orienta- 
tion are made using the transverse energy flow measured in the target calorimeter (TCal), 
and participant calorimeter (PCal). Both calorimeters have 27r azimuthal coverage and, 
combined together, they provide nearly complete polar angle coverage as viewed from the 
nucleus-nucleus center of mass system: TCal and PCal cover the pseudorapidity regions 
—0.5 < r] < 0.8 and 0.8 < r] < 4.2, respectively PJT9[]. The pseudo-rapidity r] = - \n 
tan(^/2) is defined in terms of the polar angle 6 in the laboratory frame. 



be found there. Here, the distributions d^N/ (ptdptdydcj)) of identified protons and charged 
pions are generated in a coordinate system where the x- and z-axes span the reaction plane. 
The azimuthal angle in this system is defined as = (piab — ^^i, where \E'i is the reaction plane 
angle, measured for every event using the direction of the transverse energy flow in TCal and 
PCal, and (piab is the azimuthal angle of an individual particle in the laboratory frame. The 
X-axis is defined in such a way that it points in the direction of the transverse energy flow at 
forward rapidities. The transverse momentum components Px and Py of identified particles 
are evaluated in this coordinate system, and the particle distributions with respect to the 
reaction plane are constructed. The azimuthal anisotropy of particle production is studied 
by means of Fourier analysis of azimuthal distributions [|T|pl,|T2|. This yields the rapidity, 
transverse momentum, and centrality dependence of the Fourier coefficients f„ (amplitude 
of n-th harmonic) in the decomposition: 



III. ANALYSIS 



The structure of the analysis is very similar to that used in and many details can 



E 



d^N d^N 1 d'^N 



(1 + 2v[ cos(</)) + 2v'^ cos(2<^) + 2v'^ cos(3</)) + ...). 




d^p ptdptdydcj) 27r ptdptdy 
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The E877 spectrometer provides full 271 acceptance in the spectrometer only for a very 
limited pt range (approximately pt < 50 MeV/c), but the triple differential multiplicity is 
constructed in full, using the 27r range of the reaction plane angle distribution. Note that 
the coverage of the calorimeters used for the reaction plane determination does not overlap 
with the spectrometer coverage, and thus the analysis is largely free from problems related 
to auto-correlations. 

Similarly to the analysis presented in |0, the reaction plane angle is determined 
in four non-overlapping pseudorapidity windows. The 'reaction plane resolution', i.e. the 
accuracy with which the reaction plane orientation is determined, is evaluated by studying 
the correlation between flow angles determined in different windows. Finally, the flow signals 
are corrected for this resolution. Details of this procedure are described in 0. In short, the 
true value of each Fourier coefficient Vn is obtained by: 

Vn = </(cos(n(^i - -^r)), (2) 

where v'^ is the observed signal, and (cos(?2(\E'i — "^r)) is the mean cosine of the difference 
of the reconstructed and true {"^r) reaction plane angles, characterizing the reaction 
plane resolution. The values of (cos(r;,(\E'i — ^i?)) are evaluated as outlined in As a result 
we obtain the true azimuthal distribution with respect to the reaction plane orientation: 
(fN 1 (fN 

E-— = -— (1 + 2vi cosi^iab - ^fl) + 2V2 cos(2(0;,b - ^Hr)) + ...). (3) 

a-^p In ptdptdy 

The E877 spectrometer has a relatively small azimuthal coverage. Due to this, the 
analysis of azimuthal anisotropics is rather sensitive to biases of different kinds, which could 
simulate an event anisotropy and propagate to the final results. For instance, during the 
off-line analysis, it was found that a small fraction of the recorded events (< 2-4%) have an 
anomalously high number of hits in the drift chambers. In these events the hit density was 
too high to perform reliable tracking and the high occupancy of the spectrometer was found 
to be correlated with the orientation of the reaction plane. A bias due to this was avoided 
by removing these events from the analysis completely and not only from the sample used 
to generate spectra. The tracking efficiency in the spectrometer (typically on the order of 
90%) depends slightly on the spectrometer occupancy. A special correction for this effect 
was developed and checked by a Monte-Carlo simulation. The correction is based on the 



weighting of each track in accordance with the local track densities in the drift chambers 
and time-of-flight wall. Due to different gain factors and dead towers in the calorimeters the 
distribution in the reconstructed reaction plane angle is generally not flat. Special precaution 
was taken to make the reaction plane angle distribution as flat as possible and to remove 
possible biases (see Appendix A). 

The Fourier coefficients of azimuthal distributions evaluated by the procedure described 
above (and corrected for the reaction plane resolution), are combined with the measurements 
of pt and rapidity spectra [T3|. In this way the triple differential distributions in pt, and 



(see Equation (|^)) are determined and analyzed. 



IV. RESULTS 



Directed Flow 



Anisotropic flow reveals itself already in the dependence on the azimuthal emission angle 
of the inverse slope parameter of the invariant spectra Ed^N/ d^p = d^N/ {ptdptdydcj)) = 



d^N/ {mtdmtdydcp) when plotted as a function of rrit — m for a given 0. Here, rrit = ypt + m? 
is the transverse mass. We extract the inverse slope parameter TB{<f)) by fitting the spectra 
in the region rrit — m > 0.1 GeV/c^ by a thermal (Boltzmann) function 

^d^N ( rrit — m\ 

i?^ocm,exp(^-^^j. (4) 

The shape of the spectra is not perfectly reproduced by this function, and in order to obtain 
a good quality description of the spectra in the entire region, the weights of all points 
were chosen to be equal for the fit and not in accordance with statistical errors. 

In Fig. 1^ we show the rrit spectra of protons emitted into the rapidity interval 2.8 < 
y < 2.9 at different angles relative to the reaction plane together with these thermal fits for 
different centralities. Clearly, the spectra of protons emitted in the direction of flow (0=0) 
are signiflcantly flatter than those of protons emitted in the opposite direction (0 = tt). To 
visualize the angular dependence of the inverse slope parameters the results from the flt are 
shown in Fig. ^ for three rapidity intervals. Here Tg is plotted as a function of the azimuthal 
emission angle. The results are presented for four centrality regions, selected in accordance 
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with PCal Et and corresponding to the values of atopiEx) / dgeo ~ 23-13%, 13-9%, 9-4%, 
and <4% (see Fig. 4 in 0). The value of atopiEx) is obtained by an integration of da/dET 
from a given value of Et to the maximal one observed, and the geometric cross section is 
defined as a^eo = T^rl{A^/^ + A^/^f = 6.13 b, with A = 197 and ro=1.2 fm. 

The results shown in Figs. |^ and |^ are not corrected for the reaction plane resolution. 
The correction would increase the difference between the maximal and minimal values of 
Tb by about a factor of 1.5. Such a correction (in terms of Tg) is rather complicated; for 
the quantitative description of the flow effects we use a Fourier analysis of the azimuthal 
distributions [^^|12| where all corrections are implemented. 

We quantify directed flow by Vi, the amplitude of the first harmonic in the Fourier de- 
composition of the azimuthal particle distribution defined in Equation (|1|) and corrected for 
the reaction plane resolution. The coefficient vi is analyzed as a function of transverse mo- 
mentum for different rapidity bins and collision centralities. The results for Vi{pt), corrected 
for the reaction plane resolution, are presented in Figs. for protons and charged pions. 

The error-bars shown in all figures represent statistical errors only. The systematic un- 
certainties have mostly two sources: i) The uncertainty in the determination of the reaction 
plane resolution leading to a relative error in vi of the order of 5-10%, similar for all particle 
species, ii) The uncertainty in the occupancy correction, the accuracy of which we estimate, 
by inspection of Monte-Carlo simulations, to be of the order of 20-30%. The correction 
itself is different for different particle species. It is negligible for positive pions (and for 
negative pions from the 1993 run), which are registered in the low occupancy region of the 
spectrometer. The correction is maximal, and reaches absolute values of about 0.1 for the 
data shown, for protons at low pt and/or high rapidities. The rapidly increasing uncer- 
tainty in the occupancy correction in the spectrometer region close to the beam limits our 
measurements of proton flow at very low pt- Multiplying the uncertainty in the occupancy 
correction with its absolute value we end up with a systematic error in the absolute value 
of Vi of about 0.03 where the correction is maximal. 

As can be seen from Fig. ^, proton emission is very strongly correlated with the orienta- 
tion of the reaction plane. Protons of larger pt have larger values of vi with some tendency 
to saturation in the highest pt region. The largest flow signal observed corresponds to a dif- 
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ference in the high pt proton yields along (0 = 0) and opposite to the flow direction (0 = vr) 
of almost a factor of 10 (as can also be seen in Fig. Both positive (Fig. ^ and negative 
(Fig. ^ pions exhibit weak flow in the direction opposite to that of protons (negative values 
of Vi) over most of the pt region studied. The maximum negative values of vj are about 
—0.1, signiflcantly less in magnitude than the flow signal observed for protons. [| 

Independent of any flow scenario fi(pt) must vanish at pt = 0, due to the continuity of 
the spectra. The small non-zero values of Vi for the lowest pt bin in Figs. |^ and |^ (for low 
centralities and rapidities close to beam rapidity; positive for negative pions and negative for 
positive pions) are mostly due to the flnite bin size of the data. The details of the behavior 
of f 1 at very low pt are presented in Fig. |^ where we show, in the same plot, the results for 
positive and negative pions in the rapidity region 3.2 < y < 3.6 for two different centralities 
and with much flner binning. 

Not only are the sign and magnitude of the flow signal very different for protons and 
pions but also the functional dependence on pt varies with particle species. The Vi values 
for protons grow almost linearly with pt over the entire pt region. Inspecting the transverse 
momentum dependence of Vi for pions one can distinguish three regions in pf. i) The very low 
Pt region (below pt ~ 0.1 GeV/c) where in fact the flow signals of positive and negative pions 
are different. Positive pions show vi values decreasing rapidly and monotonously towards 
negative values. Conversely, for negative pions, the flow signal becomes at flrst positive, 
reaches a peak at about pi=0.01-0.02 GeV/c, and then decreases towards negative values 
where the flow signals for both pion charges merge. This is best seen in Fig. 0. The merging 
point appears to depend on centrality, moving to lower pt for the more central events, ii) An 
intermediate region, approximately at 0.1 < < 0.3 GeV/c, where Vi is negative and only 
weakly dependent on pt for both pion charges. And flnally iii) the high pt region where f i 
begins to rise and becomes positive. There is in fact a very systematic rapidity dependence 
of the zero-crossing point: It moves to lower pt values with increasing rapidity occurring at 



Due to a "hole" in the pt acceptance of negative pions for the magnetic field polarity used during 
the 1994 run we have combined data from the 1994 run and from the lower statistics 1993 run 
(shown in Fig. |6| as open points), where data exist for both polarities of the magnetic field. 
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= 0.5-0.6, 0.4-0.5, and 0.3-0.4 GeV/c for the three rapidity bins at ?/=2.8-3.2, 3.2-3.6, and 
3.6-4.0, respectively. 

The centrahty dependence of all flow signals is rather pronounced; in the analyzed cen- 
trality region the magnitude of flow for all particles decreases for more central collisions. 
Directed flow of beam rapidity protons shows a relatively smaller centrality dependence 
than that of protons at lower rapidities (see also Fig. ^ below) . 



Higher Harmonics in Azimuthal Distributions 

Higher order harmonics {v2 and ^3) in the particle azimuthal distributions have also been 
analyzed. The accuracy in the evaluation of the contribution of higher harmonics deteriorates 
with increasing order due to the flnite reaction plane resolution (see the discussion in |]^). 
This results in larger relative errors. Our results for the proton elliptic flow (amplitude 
of the second harmonic in the proton azimuthal distributions) as a function of transverse 
momentum are presented in Fig. ^ A clear positive signal is observed in the high pt region 
for rapidities of 2.6 and larger. There is an indication that the signal moves to higher pt 
values with decreasing rapidity. This combined with the smaller acceptance in pt at rapidities 
below 2.6 may be the reason that no significant signal is observed there. In our previous 
measurements of transverse energy and charged particle flow a clear signal of elliptic flow 
was observed at all values of pseudorapidity. The observed positive values of f 2 correspond 
to an elliptically shaped distribution with the major axis lying in the reaction plane. This 
orientation of the elliptic component is perpendicular to what was measured at lower beam 



energies 



Within our spectrometer acceptance, pions do not exhibit elliptic flow at the level of 
2-3%. 

The (absolute) accuracy in measuring ^3 is about 0.1. For all particles, rapidities, and 
centralities of the collision the observed signals are consistent with zero within this accuracy. 
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Mean Directed Transverse Momentum 



Weighted with pt and its probabihty distribution the coefficient vi yields (px), the mean 
value of the transverse momentum projected onto the reaction plane: 

Like any other integral quantity, (px) contains less information than vi. We nevertheless 
calculate this quantity in order to compare our results with results from other experiments 
and model predictions. This is done by using spectra [T^ measured with the same apparatus. 
Our results of the value of (px) for protons at different centralities are shown in Fig. Due to 
the experimental acceptance in pt, a model-independent evaluation of this quantity is possible 
only at rapidities y > 2.8. Where it becomes necessary, we extrapolate dN/dpt to high pt 
using a thermal parameterization (as used e.g. in Equation (D). The filled points in Fig. ^ 
correspond to {px) calculated in accordance with Equation (|]), using the parameterization 
of Vi{pt) shown in Fig. ^ as the solid (upper) curves to extrapolate vi into the pt range not 
measured (for an analytic expression of this parameterization see Equation (|ll|) below). The 
small difference between the parameterization and the data in the low pt region does not 
contribute visibly to the final result for {px)- The contribution of the non-measured high p^ 
part of the spectra to {px) is relatively small at rapidities y > 2.8 (less than 10%), where 
the acceptance in pt is large. At rapidity y = 2.5 this contribution accounts for about 40% 
of the value of (px)- The error-bars shown do not include the uncertainty in the effective 
slope parameters used for the extrapolation of the pt spectra or any systematic uncertainty 
associated with the extrapolation of Vi. 

In order to assess the systematic error we evaluate the same quantity {{px)) by using 
a different parameterization of the invariant triple differential distribution with (px) as a 
parameter. We use the following functional form to parameterize d^N/d^p of protons in a 
given rapidity bin: 

d'N/d'p<xm[exp(-^^i^Y (6) 



m[ = ^{px - {Px)y +Pl + m^, (7) 
i. e. a thermal distribution with respect to an origin displaced along the p^^-axis. With the 

10 



effective slope parameters taken from [|T8[ we use function to fit the experimental values 
of Vi{pt) (see Fig. ^ dashed (lower) curves). The extracted values of (p^) are shown in 
Fig. P as open symbols. The difference between the results obtained with the two different 
parameterizations (filled and open symbols in Fig. ^ gives an idea of the systematic un- 
certainty of the results, not including the systematic uncertainty in vi itself in the range 
where it is measured (see above). The latter could be important for the lowest rapidity 
region of 2.2 < y < 2.4, where the flow signal is very small; there we estimate an associated 
systematic uncertainty of {px) of about 20 MeV. 

The evaluation of (px) values for pions was done by convoluting the experimental values 
of Vi{pt) with the spectra without any extrapolation to the high pt region. This is possible 
due to the relatively large pt acceptance for pions in the rapidity interval studied. The 
extracted values of (px) for pions are about an order of magnitude smaller than those for 
protons. The results for the centrality region atop/cygeo ~ 9-13% (centrality 2) are presented 
in Fig. [TT] (squares and triangles) together with the corresponding proton results and will 
be discussed below. 



V. DIRECTED FLOW DISCUSSION 

Moving thermalized source vs. absorption 

One of the simplest pictures of directed flow would be the motion in the transverse 
plane of a thermalized source, localized in rapidity. Assuming for the simplest case of no 
radial expansion, the invariant momentum spectrum for particle emission from a thermalized 
source is described by: 

1 dN 



oc E*e-^*/^, (8) 



rat dmtdydcj) 

where T is the temperature, and E* is the particle's energy in the rest frame of the source. 
In the case of a thermal source moving in x-direction with velocity (3x (which we call the 
directed flow velocity), the value of E* can be obtained by a Lorentz transformation: 

E* = -fE- (3xixPtCos{<i)), (9) 
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where 7^; = 1/y 1 — and the energy E is evaluated in the system moving longitudinally 
with the same velocity as the source, E = mtCosh.{y — y*). Here y* is the source rapid- 
ity. Using Equations (H) and (||) one can evaluate the first Fourier coefficient f 1 by direct 
integration to yield 

^ h{o-{cm-h{o)T/E 

loiO-^hiOT/E ' ^ ^ 

where ^ = PtPxIx/T, and Io{^), hiO ^^re the modified Bessel functions. For protons, rrit^ T 
and f 1 does not depend on the, generally unknown, value of y* . Since we expect that /J^, -C 1 
it is useful to note that in that case, for relatively small values of pt, vi depends almost 
linearly on pt'- 

„,«M) {^pM, (11) 

/„«) \ 2T j ^ ' 

The solid lines shown in Fig. ^ correspond to fits to the data using Equation (pTD and inverse 



slope parameters Tb from as T. Overall, the fit is rather good. Taking into account that, 
in general, one cannot interpret the inverse slope parameters Tg as a source temperature 



we would like to note that Equations (p!0| , [TTD are still valid for the case that the invariant 



spectra at a fixed rapidity have a thermal shape with some effective temperature constant. 

In Fig. the extracted values of (3x of the proton source are shown for different collision 
centralities as a function of proton rapidity. One can see that the transverse source velocity 
grows with increasing rapidity, possibly peaking around rapidity y=3.1. The maximum 
values are about 10 % of the speed of light. This transverse velocity is found to decrease 
with increasing centrality of the collision. These results have to be put into perspective with 
the transverse expansion velocities fitted to spectra of various particle species in the same 
reaction and close to mid-rapidity: in such an analysis transverse expansion velocities of 
about 50 % of the speed of light are required to describe the data |T5|. This effect is largest 



in central collisions. The transverse expansion velocity decreases away from mid-rapidity. 

Looking more closely at the fits in Fig. ^ one notices a small deviation of the fit using 
Equation ([TT|) from the data at small transverse momenta {pt < 0.1 GeV/c) and at rapidities 
y < 2.8. There, protons exhibit a weaker or even opposite fiow than expected from the simple 
model discussed above. This could be an indication of transverse expansion as pointed out in 



reference [^. There, it is shown that, depending on the relative magnitude of the sideward 
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flow velocity (3x and the transverse expansion velocity /3t, a reduction or even sign change of 
f 1 at small pt is possible. 

For pions we observe a dependence of Vi on pt very different from that for protons. It 
implies that the physics of pion flow is different from that of a moving source alone (just 
as one would expect). A possible explanation of the pion flow signal could be found in a 
superposition of different effects: Absorption in nuclear matter, Coulomb interaction with 
comoving protons, and sideward motion of the source. 

For pions produced in the center of the collision volume and moving in forward direction 
with velocities close to speed of light, one would expect more nuclear matter on the side 
characterized by the direction of nucleon flow (positive x in our notation) and, consequently, 
more absorption. The effect of pion absorption on nucleons with energy-, or for a fixed 
rapidity, pt-independent cross section would lead to a negative and constant value of vi{pt). 
Such absorption should be comparable for positive and negative pions. On the other hand, 
it is known that in the relevant momentum range the pion-nucleon elastic and total cross 



sections are strongly peaked at the A resonance In a frame where the nucleon is at 

rest the cross section peaks at a pion momentum of 0.3 GeV/c and falls off rapidly for 
larger momenta. If the absorption of pions would occur on nucleons of the same rapidity, 
this would mean that the negative values of vi should increase towards zero for pt > 0.3 
GeV/c. Most of the pions we observe are in a rapidity range where there are very few 
nucleons (forward of beam rapidity). Indeed, analysis of the proton rapidity distribution 
indicates that nucleon sources (fireballs) are distributed evenly over plus and minus one 
unit of rapidity in the cm. frame, i.e. that the most forward nucleon source is at rapidity 2.5 
in the laboratory. For pions at rapidities 3.0, 3.4, and 3.8 this would imply that absorption 
effects should decrease for transverse momenta larger than 0.30, 0.27, and 0.20 GeV/c. The 
data indeed show this trend with rapidity at about the pt values expected (see Figs. ^||, and 
0)- 

At the same time Coulomb interaction with positively charged nuclear matter would be 
different for pions of different charge, and could qualitatively explain the difference in flow 
signals at very low pt shown in Fig. |^. Negative pions are attracted to the protons leading to 
positive values of Vi (as argued in the case of absorption on nucleons above) while positive 
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pions are repelled leading to more negative values of vi. 

The rise of Vi for pions towards large pt could reflect the sideward motion of the source, 
the same way it does for protons. 



Comparison with the results on Nc and Et flow 



Fig. 11 shows, for the Ex-Tange 200-230 GeV, the mean directed transverse momenta for 
protons, positive, and negative pions. Due to the acceptance of the spectrometer the proton 
and pion results are forward of rapidity 2.2 and 2.8, respectively. Results were obtained 



using Equation ( [Tl| ) to complete the integral (||) beyond the range where Vi is measured. 



The error-bars shown include an error in the effective temperature of about 7 % in addition 
to the statistical errors. The symmetry of the collision system allows reflection of these 
results about mid-rapidity, shown in Fig. ^ by the open symbols, providing values for 
protons backward of rapidity 0.9 and for pions backward of rapidity 0.3. To flll the gap at 
intermediate rapidities we use our measurement of the flow signal in transverse energy Et 
and charged particle multiplicity Nf. as published in 0. In this paper transverse energy 
and charged particle flow were decomposed into nucleon and pion flow under a few simple 
assumptions. Both were studied as a function of pseudorapidity. Using these results and a 
simple parameterization of the proton and pion spectra in y and pt from the AGS experiment 
E866 P2[ and our own measurement |]T3 which together provide complete rapidity coverage. 



we evaluate (px) as a function of rapidity. More speciflcally, we parameterize the spectra by 
expressions of the form: 

d^N {y-yof mt-m 

oc exp( — ) ■ mt exp( ), (12) 



mtdmtdy 2 ■ Tb{y) 



y 

and 



'T 

with values of ay = 0.89 (0.88), To = 0.26 (0.15) GeV, and ar = 1.06 (1.70) for pro- 
tons (charged pions). Using these parameterizations we calculate the mean rapidity y{r)) 
(weighted with Et) and the width of the rapidity interval associated with a bin in pseudo- 
rapidity. Then, under the assumption that ff"^ and Vi* are similar, one obtains (px) as a 
function of rapidity: 
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Tb{y) = To exp( (13) 



{Px{y)) ~ {pt{y)) -vi^iy). 



(14) 



The resulting mean transverse momenta are presented as stars in Fig. together with our 
spectrometer results, for protons and pions; the errors shown are a propagation of the sta- 
tistical errors in the measured flow of Et and N^. The horizontal error-bars shown with the 
stars represent the widths of the rapidity region actually contributing to the pseudorapidity 
bins for which Vi was measured. 

Overall, the agreement between the results of the current analysis of spectrometer data 
at forward (and backward) rapidities and the results derived from the measurements of Et 
and Nc flow is good. For the latter method we had estimated the systematic error for proton 
flow to be about 10% 0; this can also be judged from the symmetry of the data shown in 
Fig. |TT] with respect to midrapidity (note that all stars correspond to measured points). The 
dotted line has been added to aid inspection of the data for the necessary antisymmetry. 
The first and the next to last proton point have equal distance from mid-rapidity and should 
be identical. They differ by 10 MeV/c (out of 80 MeV/c). For the spectrometer results the 
systematic errors are largest at the lowest rapidities shown (see Fig. ^ and discussion) where 
they are estimated to be 20 MeV/c. 

Within these errors the points obtained from the two methods merge smoothly and taken 
together they provide a measurement of (p^) for protons and charged pions as a function 
of rapidity over practically the entire relevant rapidity range. One can see a proton flow 
signal rising away from midrapidity to maximal values of about 130 MeV/c close to beam 
and target rapidities. The pion signal is in comparison very small, about 1/20 or the proton 
signal, but nevertheless significant. It is directed opposite to the nucleon flow signal. 



Comparison with RQMD predictions 



In Fig. |TT] we compare our results for (px) with the predictions of the RQMD event 
generator p!7|JIO|] , versions 1.08 and 2.3. Version 1.08, run in cascade mode, has been used for 



comparison to our earlier results for Au + Au collisions on pseudorapidity distributions p3 



proton and pion rapidity [0 and transverse mass [jT5[ distributions as well as results for 
Si+A and p+A collisions. For the Au + Au system it was found that overall many features 
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of the data are predicted correctly with one striking deviation: generally proton spectra are 
significantly too steep in the model; at midrapidity the inverse slope constant is only 2/3 of 
the experimental value. Predictions from this version of the model have also been compared 
to our experimental results for fiow in Et and Nc and generally the model exhibited 
too little fiow as compared to the data. Comparing to the proton data in Fig. |TT], we find 
that the model also underpredicts the proton fiow by nearly a factor of two. For pions, the 
sign change of the fiow relative to protons is properly predicted but in magnitude the model 
overpredicts this opposite pion fiow by a factor of 2 to 3. 

Meanwhile, a new version of RQMD (version 2.3) has been developed ||10|- It has been 
known for some time that taking into account mean field effects by simulating in the 
model a Skyrme-type nucleon-nucleon potential increases the slopes of proton spectra and at 
the same time increases the proton fiow signal while reducing the opposite pion fiow signal 
in magnitude. We have used the most recent version of the model in the so-called "mean 
field" mode to compare to the present data as well. Indeed, in this version of the code the 
proton fiow signal is larger by about 50% for the proton (p^) at the peak and the data are 
described rather well (see Fig. |1^). Simultaneously, the opposite pion fiow is reduced, again 
to a level consistent with the experimental data. It should be noted that at the same time 
the slope difference in the proton spectra is reduced (at midrapidity the slope constant of the 
model is now 3/4 of what is observed in the data) but a deviation persists at all rapidities. 

Fig. ^ shows that the agreement of the RQMD predictions (version 2.3, mean field) 
with the experimental results for (p^.) of protons as displayed in Fig. |lT] may indeed be 
accidental. The functional dependence of the proton fiow signal, expressed now in terms 
of the first Fourier coefficient Vi as a function of pt, is very different from what is seen in 
the data. In the model, the flow signal rises at flrst with increasing pt and then becomes 
rather flat. In fact, in the rapidity bin around y = 2.9 the model overpredicts vi but at the 
same time underpredicts the overall slope of the proton spectrum by 25 % giving overall 
agreement in (px)- An interesting question arises whether the rapidly rising and saturating 
Pt dependence as exhibited by the RQMD predictions (which is in fact a general trend 
observed for all rapidity bins) is due to the rescattering-type mechanism which produces 
flow in a cascade code in contrast to a hydrodynamic flow mechanism (reflected e.g. by the 
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Pt dependence of Equations (|TO|,|TT])). 

The elliptic flow signal seen in the anisotropy of Et and Nc was found to be in good 
agreement with the RQMD (version 1.08) prediction in our previous publication [0. As 
discussed above, the sign of this elliptic flow signal changes from negative (preferential 
emission perpendicular to the reaction plane) at lower beam energies (1-2 GeV/nucleon 
range at the Bevalac and SIS) to a positive signal (preferential in-plane emission) at the 
present energy. It was recently noted by Sorge |p that the V2 anisotropics are sensitive to 
the pressure at the maximum compression. Positive values of f 2 were in fact predicted in 



a hydrodynamic model by Ollitrault for very high energies []lT|. The observed change in 
sign between Bevalac/SIS and AGS energies therefore implies that at AGS energies already 
the final in-plane flow overwhelms the initial shadowing effect. It would be interesting to 
deduce, from this information, the pressure achieved at maximum baryon density in Au + 
Au collisions at the AGS. 



VI. SUMMARY AND CONCLUSION 



The measured flow signals of identified particles show that, in semi-central collisions, 
protons (nucleons) exhibit strong directed flow. Pions of low pt exhibit flow in the opposite 
direction; at higher pt pions start to flow in the same direction as protons. The results 
on (px) values for proton and pion are complementary in acceptance and match well with 
nucleon and pion flow values derived from the same experiment by measurement of Ej- and 
Nc flow 10. Taken together, these results represent a directed flow measurement of protons 
and pions over nearly the entire rapidity region in Au+Au collisions at 11 A GeV/c. The flow 
is found to be maximal around beam and target rapidities with values of (p^) for protons of 
about 130 MeV/c and 5% of that for charged pions. 

The observed positive values of the second harmonic amplitudes of the proton azimuthal 
distributions correspond to preferential particle emission in the reaction plane. This obser- 
vation agrees with our previous measurements of E^ and A*";, elliptic flow f^. 

The nearly linear dependence of the proton flow signal vi on pt can be interpreted in the 
framework of a transversely moving thermal source. The corresponding source velocity jS^ 
appears to reach values of 0.1 in the beam rapidity region and for semi-central collisions. 
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The more complicated pt dependence of the pion flow signals indicates that the effect there 
is probably a superposition of several effects such as absorption, Coulomb interaction and 
overall motion of the source (as for protons). 
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APPENDIX 



Flattening of the Reaction Plane Distribution 

We apply the following procedure to correct for the non-flatness of the reaction plane 
distribution. Any fiattening of the distribution means a correction to the reaction plane 
angle. We explicitly introduce this correction defining a new angle as: 

% = ^1 + A^i. (15) 

where A^i is written in the form: 

A*i = ^(^„cos(n*i) + S„sin(n*i)) (16) 

n 

Requiring the vanishing of the n-th Fourier moment of the new distribution, the coefficients 
An and Bn can be evaluated by the original distribution. 

S„ = -(cos(n*i)), (17) 

Th 

An = --(sin(n*i)), (18) 

Th 

where the brackets refer to an average over events. This gives: 
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= + ^ -( -(sm(n^'i)) cos(n^'i) + (cos(n^'i)) sm(n^'i)) (19) 

n ^ 

In practice, we flatten the reaction plane distribution up to the fourth harmonic (n=4). 
Note that, due to the small values of and (typically of the order of a few percent), 
such a flattening of the distribution docs not have any effect on the reaction plane resolution. 
It can also be shown that the same flattening procedure removes possible trigger biases (due 
to imperfect calibration, dead channels or any other asymmetry) at least up to the second 
order. 
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FIGURE CAPTIONS 



1. The E877 apparatus. 

2. Proton nit distributions in the rapidity interval 2.8 < y < 2.9 for different centrahties 
(indicated by an Et range in GeV) together with a fit using the function given in 
Equation (^). The inverse slope values and their statistical errors are shown in GeV. 
The open (filled) squares correspond to — 7r/4 < < 7r/4 (37r/4 < < 57r/4). 

3. The dependence of the inverse (Boltzmann) slope Tg of the proton distributions as a 
function of azimuthal angle. The open points are reflections of the filled points about 
= 0. 

4. The transverse momentum dependence of the first moment (vi) of the proton azimuthal 
distributions for different particle rapidities and centralities of the collision. The solid 
and dashed curves are fits using functions given in Eqs. (|TTp and (^, respectively. For 
a description see Sections [I^ and 0. 



5. Same as Fig. ^for tt"*". 

6. Same as Fig. ^ for tt~ . The open and filled symbols correspond to data from the 1993 
and 1994 runs. 

7. Fourier coefficients Vi for pions in the low pt region for the rapidity interval 3.2 < 
y < 3.6 and two different centralities. Solid symbols represent data for positive pions. 
Open symbols are for negative pions (circles and squares correspond to results from 
the 1993 and 1994 runs, respectively). 

8. The transverse momentum dependence of the second moment (^2) of the proton az- 
imuthal distributions for different particle rapidities and centralities of the collision. 

9. The mean projection (p^) of the proton transverse momentum onto the reaction plane 
as a function of rapidity for different centralities of the collision. Results for the second 
and fourth centrality bin are shifted to the left by 0.025 for clarity of the picture. Solid 
and open symbols correspond to fits of Vi{pt) with functions given in Eqs. (0) and (^, 
respectively. 

20 



10. The transverse velocity (3x of a thermal (proton) source for different rapidities and 
centralities of the collision. For details see text. 

11. Mean projection {p^) of the proton and pion transverse momentum onto the reaction 
plane as a function of rapidity for the centrality bin with Et = 200-230 GeV. Solid 
circles, squares, and triangles correspond to measurements in the spectrometer. Open 
circles, squares, and triangles are reflections of the filled symbols about midrapidity. 
The results derived from measurements of Et and Nc flow for nucleons and pions 
are shown by stars (see text). The dotted line is added to aid inspection of the solid 
stars for antisymmetry about mid-rapidity. Results from calculations using 2 versions 
of RQMD (1.08 in cascade mode, 2.3 in mean field mode) are depicted as histograms. 

12. Comparison of measured values of vi for protons in two rapidity bins and the centrality 
bin with Et = 200-230 GeV with those predicted by 2 versions of the event generator 
RQMD. 
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FIG. 1. The E877 apparatus. 
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FIG. 2. Proton rrit distributions in the rapidity interval 2.8 < y < 2.9 for different centralities 
(indicated by an Et range in GeV) together with a fit using the function given in Equation (Q). 
The inverse slope values and their statistical errors are shown in GeV. The open (filled) squares 
correspond to — 7r/4 < < 7r/4 {Sir/A < (j) < bir/A). 
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FIG. 3. The dependence of the inverse (Boltzmann) slope Tb of the proton distributions as a 
function of azimuthal angle. The open points are reflections of the filled points about ^ = 0. 
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FIG. 4. The transverse momentum dependence of the first moment (vi) of the proton azimuthal 
distributions for different particle rapidities and centralities of the colhsion. The sohd and dashed 
curves are fits using functions given in Eqs. (|ll|) and (|6|), respectively. For a description see 
Sections IV and |V[ 
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FIG. 5. Same as Fig. || for tt^ 
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FIG. 6. Same as Fig. |^ for vr . The open and filled symbols correspond to data from the 1993 
and 1994 runs. 
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FIG. 7. Fourier coefficients vi for pions in the low pt region for the rapidity interval 3.2 < y < 3.6 
and two different centrahties. Sohd symbols represent data for positive pions. Open symbols are for 
negative pions (circles and squares correspond to results from the 1993 and 1994 runs, respectively). 
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FIG. 8. The transverse momentum dependence of the second moment {V2) of the proton 
azimuthal distributions for different particle rapidities and centraUties of the colhsion. 
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FIG. 9. The mean projection {p^) of the proton transverse momentum onto the reaction plane 
as a function of rapidity for different centrahties of the cohision. Results for the second and fourth 
centrality bin are shifted to the left by 0.025 for clarity of the picture. Solid and open symbols 
correspond to fits of vi{pt) with functions given in Eqs. (lll|) and @, respectively. 
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FIG. 10. The transverse velocity Px of a thermal (proton) source for different rapidities and 
centralities of the collision. For details see text. 
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FIG. 11. Mean projection {px) of the proton and pion transverse momentum onto the reaction 
plane as a function of rapidity for the centrahty bin with Et = 200-230 GeV. Sohd circles, squares, 
and triangles correspond to measurements in the spectrometer. Open circles, squares, and triangles 
are reflections of the filled symbols about midrapidity. The results derived from measurements of 
Et and flow |Q| for nucleons and pions are shown by stars (see text). The dotted line is 
added to aid inspection of the solid stars for antisymmetry about mid-rapidity. Results from 
calculations using 2 versions of RQMD (1.08 in cascade mode, 2.3 in mean field mode) are depicted 
as histograms. 
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FIG. 12. Comparison of measured values of vi for protons in two rapidity bins and the centrality 
bin with Et = 200-230 GeV with those predicted by 2 versions of the event generator RQMD. 
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